Dynamics of photoinduced Charge Density Wave-metal phase transition in K0.3M0O3 
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We present first systematic studies of the photoinduced phase transition from the ground charge 
density wave (CDW) state to the normal metallic (M) state in the prototype quasi- ID CDW system 
K0.3M0O3. Ultrafast non-thermal CDW melting is achieved at the absorbed energy density that 
corresponds to the electronic energy difference between the metallic and CDW states. The results 
imply that on the sub-picosecond timescale when melting and subsequent initial recovery of the 
electronic order takes place the lattice remains unperturbed. 



In superconductors (SC), it has been known for 
decades that an intense laser pulse can non-thermally de- 
stroy the superconducting ground state The energy 
required to destroy superconductivity should, in case all 
the absorbed optical energy is kept in the electronic sub- 
system during the process of superconductivity suppres- 
sion, be equal to the condensation energy (energy differ- 
ence between the free energy of the SC and normal states 
at T = K) . In recent experiments on MgB2 Q and the 
high temperature superconductor L^-zSr^C^C^ Q it 
has been shown, however, that the absorbed optical en- 
ergy required to suppress superconductivity is substan- 
tially higher than the thermodynamically measured con- 
densation energy |3[. This discrepancy was accounted 
for by considering in detail all energy relaxation path- 
ways on the timescale when superconductivity suppres- 
sion is achieved. From this analysis it follows that on 
this timescale a quasi-equilibrium between the density 
of quasiparticles and high frequency phonons is achieved 
with most of the absorbed energy density being stored in 
the phonon subsystem 0, 0| ■ 

Charge density wave systems present another broken 
symmetry ground state. Here upon cooling through the 
CDW transition temperature the translational symmetry 
is broken [J]. The appearance of the long-range charge 
density modulation is accompanied by the appearance of 
the gap in the single particle excitation spectrum at the 
Fermi level, while the collective excitations of the CDW 
state are the so called amplitude and phase mode 
While real-time studies of photoexcited quasiparticle and 
collective mode dynamics in CDW compounds have been 
quite extensive in the weak and moderate perturbation 
regime [1, [E S H @| , and a reasonable understanding of 
the underlying relaxation processes has been obtained, 
systematic studies in the high perturbation regime, where 
the energy of the optical excitation pulse is enough to 
drive the phase transition from the CDW ground state 
to the normal metallic state 1(| 11 1 are still lacking. 



In this Letter we report on the first systematic study 
of carrier and collective mode dynamics in a prototype 
quasi-1 dimensional CDW K0.3M0O3. Systematic tem- 
perature and excitation density dependent measurements 



of the photoinduced reflectivity changes reveal that the 
phase transition from the ground CDW state to the nor- 
mal metallic state can be achieved on the femtosecond 
timescale. From the energy conservation law it follows 
that the phase transition is non-thermal in origin; i.e. 
the phase transition is not a result of a simple heating 
of the sample to above the critical temperature. The 
absorbed energy density required to optically induce the 
phase transition is found to be comparable to the elec- 
tronic energy difference upon CDW condensation. These 
results give new insight in the ultrafast processes gov- 
erning the relaxation dynamics in low dimensional CDW 
systems. In particular, the results suggest that on the 
timescale shorter than the period of the characteristic 
lattice vibrations (« 0.6 ps in K0.3M0O3, which is the 
inverse frequency of the amplitude mode 
charge density modulation is suppressed while the lat- 
tice remains unperturbed keeping the 2 modulation. 

In the experiments described here, we utilized a degen- 
erate optical pump-probe technique to study the excita- 
tion intensity and temperature dependence of the pho- 
toinduced reflectivity dynamics in single crystals of blue 
bronze K0.3M0O3. We used a commercial TkSapphirc 
amplifier producing 6 fiJ, 50 fs laser pulses at A = 800 
nm (photon energy of 1.55 eV) at a variable repetition 
rate between 9 and 250 kHz. The laser was used as a 
source of both excitation and probe pulses. Samples were 
mounted in an optical helium flow cryostat, with both 
excitation and probe beam entering the sample at near 
normal incidence. Due to the strong anisotropy of the 
induced changes in reflectivity with respect to light po- 
larization (BJ, the probe laser beam was polarized along 
the chain [010] direction, while the excitation beam was 
polarized along the perpendicular [102] direction To 
ensure a homogeneous excitation profile, the diameter of 
the pump beam at the sample position was twice the 
diameter of the probe beam. To determine the photoex- 
citation density at the position of the sample with high 
precision we used a beam profiler. We fitted the beam 
profile with a Gaussian, and the excitation fluences, F, 
used throughout the paper correspond to the maximum 
fluence at the center of the beam. Low thermal con- 
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ductivity in K0.3M0O3 [14| can lead to a pronounced in- 
crease of the equilibrium temperature in the probed vol- 
ume, proportional to the average laser power. Therefore, 
the high excitation experiments were performed at a low 
repetition rate (10-30 kHz), where at F = 1 mJ/cm 2 the 
temperature increase is < 10 K. 
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FIG. 1: Photoinduced reflectivity changes in K0.3M0O3 at 
4 K following photoexcitation with 50 fs optical pulses at 
different excitation fluences [l5l ]. Inset: the F— dependence 
of the initial decay time Ti determined by fitting the data as 
described in Ref.[lq|. ti(F) displays critical behavior (dashed 
line is a guide to the eye) near the threshhold fluence of ~ 200 
(xJ/cm 2 . 

Figure 1 presents the induced reflectivity transients 
taken at 4 K at several excitation densities. At low exci- 
tation densities the data show the same behavior as pre- 
viously reported [B| . The decay dynamics of the (incoher- 
ent) electronic response show a bi-exponential decay with 
timescales t\ ~ 0.3 ps and T2 ~ 7 ps. The former one, 
showing critical slowing down upon approaching T^P, 
was attributed to the recovery of the CDW gap, while 
the second one was initially tentatively attributed to an 
overdamped phase mode [5J- Recent detailed studies of 
the dynamics as a function of F and applied external elec- 
tric field however suggest that this longer timescale more 
likely presents the second stage of the CDW recovery [lj| . 
On top of the incoherent transient an oscillatory (coher- 
ent) signal is observed whose Fourier transform, obtained 
by Fast Fourier Transform (FFT) analysis, shows several 
frequency components which can be attributed to the co- 
herently excited amplitude mode (the strongest mode at 
1.68 THz) and several other phonon modes [5, 1(1 IB]. 

While the photoinduced transient is linear in F over 
several orders of magnitude , we observe pronounced 
changes upon increasing the excitation intensity into the 
100 /iJ/cm 2 range. The electronic component shows clear 
saturation at F > 200 /iJ /cm 2 with a maximum induced 



change in reflectivity approaching 18 % [17[. One is 
tempted to ascribe this saturation behavior in reflectiv- 
ity change to the photoinduced suppression of the CDW 
order, where this difference in reflectivity corresponds to 
the change in reflectivity between the normal metallic 
and the CDW state. Since no measurement of the tem- 
perature dependence of reflectivity at optical frequencies 
is reported to date, we have performed thermomodula- 
tion measurements to determine the magnitude and sign 
of change in equilibrium reflectivity upon increasing the 
temperature to above Tj? D . We measured the reflectiv- 
ity difference between CDW and metallic state at 1.55 
eV (800 nm) by heating the sample to above the phase 
transition using a CW laser. The data show that the re- 
flectivity at 800 nm indeed decreases upon increasing the 
temperature. The change in reflectivity of ~ 10 % was 
observed upon heating from 160 K to just above T\ D . 
At temperatures above T^P the corresponding change in 
equilibrium reflectivity was less than 1 %. Therefore the 
observed reflectivity change and its saturation behavior 
upon increasing the excitation fluence arc consistent with 
the photoinduced CDW-M phase transition. 
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FIG. 2: independence of the Fourier transform of the oscil- 
latory signal at 4 K. The data were, prior to FFT analysis 
(time window of 50 ps), normalized to F for clarity. The area 
between 1.1 and 9 /xJ/cm 2 is obtained by interpolation be- 
tween the low F data recorded with the Ti:sapphire oscillator 
and the data recorded with the amplified system. Color cod- 
ing represents the FFT amplitude. Strong suppression of the 
AM (w 1.68 THz) is observed at F « 200 iiJ/cm 2 , while the 
other modes vanish at fluences that are about one order of 
magnitude higher. 

From the rise-time of the reflectivity transient it 
also follows that this transition happens on the 100 fs 
timescale after photoexcitation (rise-time is becoming 
shorter upon increasing F). The initial decay time, t\, 
shows a pronounced increase near the threshold fluence, 
followed by a rapid drop as shown in inset to Fig. 1 
(the secondary decay time ti shows only a slight decrease 
upon increasing F). This critical slowing down of relax- 
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FIG. 3: The time evolution of the FFT spectrum at 4K 
recorded at F = 12 /ij/cm 2 , 130 /xj/cm 2 , 0.38 mj/cm 2 , 
and 2.5 mj/cm 2 (panels a-d). The plots are obtained by 
time-windowed FFT analysis with a 5 ps time window. The 
AM is strongly suppressed at F > 200 /iJ/cm 2 , while the high 
frequency modes are present up to F ~ 3 mj/cm 2 . 



ation below the threshold for the CDW melting shows 
similar behavior as in the weak perturbation limit upon 
increasing the temperature towards T^P [f|. 

Additional support to the assignment of the saturation 
behavior to the photoinduced CDW-M transition comes 
from the study of the oscillatory response, shown in Fig. 
2. Here clear suppression of the amplitude mode, which 
presents a fingerprint of the CDW state, is observed at 
comparable fiuences. Figure 2 shows the 2 dimensional 
surface plot of the FFT spectrum of the oscillatory sig- 
nal as a function of excitation intensity over several or- 
ders of magnitude. Several sharp lines are observed, the 
strongest being that of the amplitude mode (AM) at 
about 1.68 THz (56 cm -1 ). The two second strongest 
modes at 2.25 THz (74 cm" 1 ) and 2.55 THz (85 cm" 1 ) 



correspond to zone- folding modes [18|, [19| . In addition, 
some weaker modes in the r ang e between 3.5-4.5 THz are 
also observed as in Raman We should also mention 
the weak side modes observed in the vicinity of the AM 
as well as the two zone-folding modes. These are not 
artefacts due to FFT. From the ratio of their amplitudes 
with respect to the main peaks, and from the consider- 
ably weaker damping constants compared to the main 
modes we ascribe these to their corresponding surface 
modes. 

While Fig. 2 clearly shows strong suppression and in- 
creased damping of the AM above F « 200 /xJ/cm 2 , sup- 
pression of other phonons follow at the excitation densi- 
ties that are about one order of magnitude higher in F. 



More information about the nature of the photoinduced 
CDW-M phase transition can be gained by looking at the 
evolution of the FFT spectrum with time at different F, 
which is presented in Fig. 3 for four different flucnccs. 
Here one clearly sees that the AM vanishes at F > 300 
/i J /cm 2 while the two zone- folding modes persist up to F 
> 2.5 mJ/cm 2 finally disappearing above F > 3 mj/cm 2 . 
The fact that high frequency modes survive above the 
threshold fluence for the photoinduced CDW-M phase 
transition suggests that on a short timescale time follow- 
ing the photoinduced CDW-M transition and subsequent 
recovery on the sub-ps timescale (inset to Fig. 1) the lat- 
tice remains largely unperturbed. In other words, pho- 
toexcitation with a 50 fs optical pulse of F > F sat induces 
melting of the electronic density modulation, which par- 
tially recovers on the sub-ps timescale, while the lattice is 
- on this timescale - uncoupled from the electron subsys- 
tem and retains its 2kj? modulation. Only in the second 
step of relaxation, which proceeds on the 10 ps timescale, 
the CDW can be described with a single order parameter 
where electrons adiabatically follow the lattice. 
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FIG. 4: a) Excitation fluence dependence of the induced re- 
flectivity maximum at several temperatures below Tj? D . The 
solid lines are fits with the simple saturation model 0]. b) 
Saturation fiuences extracted from a) and the corresponding 
absorbed energy densities in meV per unit cell volume. The 
dashed line represents the calculated T-dependence of E sa t 
given by Eq.(l), where BCS T-dependence of A was used. 

We studied the F-dependence of the amplitude of the 
electronic signal at different temperatures below T^ D 
(Fig. 4). Upon increasing the temperature saturation 
appears at decreasing values of F. In order to determine 
the temperature dependence of the saturation fluence, 
F sa t, we fit the data with the simple saturation model, 
which is described in detail in Ref. 0] . The temperature 
dependence of F sat , together with the corresponding ab- 
sorbed energy density per unit cell volume [20(, E sat , is 
presented in Fig. 4b. The first question is whether melt- 
ing of the CDW state is a result of simple pulsed heating 
of the excited volume or is it non-thermal in origin? We 
calculated the expected temperature rise AT which cor- 
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responds to E sa t, using E sat = J^° +AT c p (T)dT. Here c p 
is the total specific heat pj] ■ At T = 4 K and E sat w 60 
meV/u.c.v. we obtain AT < 40 K, while at higher To 
this value is considerably smaller (AT < 3 K at 172 K 
and E sat ks 20 meV/u.c.v.). In fact, the absorbed energy 
density required to heat K0.3M0O3 from 4K to its phase 
transition T 3 C D = 183K is about 600 meV/u.c.v., an or- 
der of magnitude higher than E sat . It follows that the 
photoinduced CDW-M phase transition is non-thermally 
driven. 

From the observation of coherently excited zone-folded 
modes at fluences up to one order of magnitude higher 
than F sat it follows, that on the sub-ps timescale after 
photoexcitation, the electrons are nearly uncoupled from 
the lattice. In this case E sat should be compared to the 
energy gain of the electronic subsystem upon CDW con- 
densation, E e i. To estimate E e i we used the mean-field 
expression in the weak coupling limit, given by Eq. 3.40 
in Ref. Q: 



(1) 



Here ti(€f) is the normal state density of states at the 
Fermi energy, ep, and A is the value of the CDW gap. 
Using A = 60 meV, ep_= 0.24 - 0.39 eV and n(e F ) 
= 4-6 eV _1 /u.c.v. H we obtain E el (4K) = 37-66 
meV/u.c.v.. This value is in excellent agreement with 
E sat (4K) as 60 meV/u.c.v., giving further support to the 
argument that during photoinduced CDW melting the 
electronic order parameter is decoupled from the lattice 
on the sub-ps timescale. At F > 3 mJ/cm 2 , which cor- 
responds to the absorbed energy density of 1 eV/u.c.v., 
other modes are also completely suppressed. This en- 
ergy density is in good agreement with the calculated 
energy density required to heat up the excited volume to 
above T 3 C D . In this regime, recovery proceeds on a much 
longer timescale which is determined by heat diffusion 
out of the excited volume. The temperature dependence 
of E sat , shown in Fig. 4b, shows good agreement with 
the expected T-dependence of E e i. The dashed line in 
Fig. 4b shows E e i(T) calculated from Eq.(l), where the 
BCS T-dependence of A was used 

We have shown that in K0.3M0O3 the photoinduced 
CDW-M phase transition is non-thermal and takes place 
on the 100 fs timescale. The good agreement between 
measured E sat and calculated E e i , the observation of the 
order parameter recovery on the sub-ps timescale, and 
the observation of zone-folded phonons high above E sa t 
suggest that during the process of melting and sub-ps re- 
covery of the electronic modulation the lattice remains 
nearly frozen. This has an important implication for un- 
derstanding the ultrafast relaxation processes in systems 
with reduced dimensionality, in particular for the sys- 
tems with strong electron-phonon interactions that lead 
to phenomena like charge density modulation. The ini- 
tial reconstruction of the CDW state is found in all sys- 



tems studied thus far to proceed on the sub-ps timescale 
S 0, H ©, EU El- Importantly, this timescale is one 
to two orders of magnitude faster than in the high-T c 
superconductors and is indeed close to the typical 
timescale for electron-phonon thermalization. The for- 
mation of the CDW requires freezing of a phonon and 
our results do imply that the lattice remains frozen in 
its modulated state on the sub-ps timescale after pertur- 
bation. Therefore, the extremely fast order parameter 
recovery in this entire class of low-dimensional materi- 
als @, S, 0, H ©, M, El could be a consequence of the 
fact that on the short timescale after photoexcitation the 
lattice remains in its unperturbed state. Thereby, the re- 
taining 2kp modulation presents a strong potential well 
driving ultrafast reformation of the charge density mod- 
ulation. Clearly, further theoretical studies as well as 
studies of the ultrafast structural dynamics are required 
to shed additional light on these fascinating phenomena. 
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